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The H.E.S.S. experiment entered its phase II with the addition of a new, large telescope named 
CT 5 that was added to the centre of the existing array of four smaller telescopes. The new 
telescope is able to detect fainter air showers due to its larger mirror area, thereby lowering the 
energy threshold of the array from a few hundred GeV down to i^(50GeV). Due to the power-law 
decrease of typical 7 -ray and cosmic-ray spectra of astrophysical sources a majority of detected 
air showers are of low energies, thus being detected by CT 5 only, which motivates the need for a 
reconstruction algorithm based on information from CT 5 alone. By exploiting such monoscopic 
events the H.E.S.S. experiment in phase II becomes sensitive in an energy range not covered by 
H.E.S.S. I and in which the Eermi EAT runs out of statistics. Eurthermore the chance of detecting 
transient phenomena like 7 -ray bursts is increased signihcantly due to the large effective area of 
CT 5 at low energies. 

In this contribution a newly developed reconstruction algorithm for monoscopic events based 
on neural networks is presented. This algorithm uses multilayer perceptrons to reconstruct the 
direction and energy of the particle initiating the air shower and also to discriminate between 
gamma rays and hadrons. The performance of this algorithm is evaluated and compared to other 
existing reconstruction algorithms. Eurthermore results of hrst applications of the algorithm to 
measured data are shown. 
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1. Introduction 

Ground-based very high energy (VHE, ~ 30GeV - 100 TeV) y-ray experiments have estab¬ 
lished a new astronomical field over the last 20 years, with more than 100 sources discovered 
already [1]. Current imaging atmospheric Cherenkov telescope (lACT) experiments like MAGIC, 
VERITAS and H.E.S.S. are able to detect the faint Cherenkov light emitted in extended air showers 
induced by VHE particles interacting with nuclei in the Earth’s atmosphere. The H.E.S.S. experi¬ 
ment [2] entered its phase II in 2012 with the addition of a hfth telescope named CT 5 to the centre 
of the existing array of four 12 m telescopes. This new telescope is the largest of its kind, with a 
mirror area of 600m^, able to detect air showers from low-energy particles. Therefore the lower 
end of the observable energy range can be decreased to ^(50GeV), compared to a few hundred 
GeV in case of the H.E.S.S. I array. Due to the power-law decrease of the number of incident par¬ 
ticles as a function of their energy a majority of all events will be detected by CT 5 only. A newly 
developed reconstruction chain called MonoReco, implementing a full event reconstruction using 
information from CT 5 only, is presented here. It is the first monoscopic reconstruction algorithm 
based on multivariate analysis techniques used within the H.E.S.S. experiment. 

2. MonoReco - A Monoscopic Reconstruction Algorithm 

Several classes of astrophysical objects are able to accelerate charged particles up to VHEs. 
Those particles, mainly protons, heavier nuclei, electrons and positrons, are deflected by inter¬ 
stellar and intergalactic magnetic helds. Thus, in the energy range accessible to lACTs, they are 
isotropised on their way to Earth. To be able to study the acceleration processes occurring in indi¬ 
vidual sources 7 rays produced by those charged particles through e.g. inverse Compton scattering 
or ;r*’ decay can be analysed. Eor a typical source the number of 7 rays detected by lACTs is lower 
than the number of detected hadrons by a factor of ^(1000). Therefore it is important to effi¬ 
ciently distinguish between those particle classes as well as to reconstruct the particles’ energies 
and incident directions. 

Most reconstruction concepts depend heavily on stereoscopy, i.e. information from at least 
two telescopes, to be able to reconstruct all of the main properties of the primary particles like 
their direction and energy as well as the type of particle. This is no longer possible for low-energy 
showers that are only seen with CT5. In order to extend the energy range accessible to H.E.S.S. 
towards the regime in which the Eermi Earge Area Telescope (EAT) [3] operates (< 100 GeV) it 
is necessary to analyse monoscopic events. Such a low-energy analysis is especially needed in the 
case of transient events like 7 -ray bursts or flaring activify from acfive galacfic nuclei because of 
fhe ~ 10"^ fimes larger effechve area of CT 5 with respect to the effective area of the Eermi EAT. 

The reconstruction chain presented here is capable of performing a full event reconstruction. 
It makes use of the so-called Hillas parameters [4], the moments of the distribution of the measured 
intensities in the individual pixels of the camera of CT 5, and multivariate analysis techniques to 
determine all relevant properties of the incident particle. 

Before multivariate analysis techniques are applied, several preprocessing steps are performed. 
Noise is removed from the images by requiring two neighbouring pixels to contain minimal inten¬ 
sities of 5 and lOphoto electrons (p.e.) or vice versa. After this image cleaning has been performed. 
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Cut parameter 

std 

loose 

extraloose 

Number of pixels 

>3 

>3 

>3 

Intensity 

> 60p.e. 

> 60p.e. 

> 35p.e. 

Nominal distance 

< 1.15° 

< 1.15° 

< 1.15° 

c 

>0.9 

>0.78 

>0.55 

0 

0.13° 

0.14° 

0.23° 


Table 1: Cut values used in the three different cut configurations introduced in the main text. 


several cuts are applied to guarantee a sufficient intensity and number of pixels. Also the nominal 
distance, i.e. the angular distance of the center of gravity (c.o.g.) of the Hillas ellipse from the 
centre of the camera, is cut on to make sure the image is not truncated at the camera edges. The cut 
values are summarised in Tab. (1) for several cut configurations. 

In each of the three main analysis steps neural networks are used to determine properties of 
the primary particles. The implementation of multilayer perceptrons (MLPs) within the TMVA 
framework [5] was used in this analysis. Before such neural networks can be used in an analysis 
they have to be trained. In the reconstruction chain presented here y-ray Monte Carlo simulations 
were used in the training stage. 

In general the shape of a y-ray air shower image in the focal plane of the telescope mirror 
is elliptical. For geometrical reasons the direction of the initial particle is a point assumed to be 
located along the major axis of the ellipse which is parameterised as a Hillas ellipse. Based on 
this assumption the direction is reconstructed using the c.o.g. together with a displacement value 5 
calculated as the distance between the c.o.g. and the direction to be reconstructed. When analysing 
data 5 is retrieved from the previously trained MLP networks working in regression mode based 
on the following input parameters: the width and length of the Hillas ellipse, the total intensity 
divided by the area of the Hillas ellipse (called density), the logarithm of the image amplitude, the 
skewness and the kurtosis of the intensity distribution in the camera. As only the absolute value of 
5 can be retrieved from the MLP network it is unclear in which direction along the major axis of 
the ellipse the 5 value has to be applied. This degeneracy can be resolved by using the asymmetry 
of the intensity distribution in the camera parameterised by the skewness. 

The reconstruction of the energy of the initial y ray is also based on MLP networks working 
in regression mode using the same input variables as the MLP networks mentioned above. 

Finally the type of the primary particle has to be determined. In order to accomplish this task 
an MLP network operating in classification mode is used. The MLP network is provided with a 
set of variables (the width and length of the Hillas ellipse, the density, the skewness and kurtosis 
as well as the length over the logarithm of the intensity) from Monte Carlo simulations in case 
of y rays and with parameters taken from so-called OFF runs, i.e. observations of regions of the 
sky containing no known sources of y rays, in case of background events. Based on such input 
events of known signal or background class the MLP network is able to construct a function which 
calculates an output parameter based on the set of input variables. This value is used to perform 
a cut in order to discriminate between yrays and hadrons. The definition of the cut value is subject 
to an optimisation procedure. Several science cases differing in the assumed spectral index of y-ray 
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differential energy spectra were accounted for by the definition of several cut sets. A std cut set 
is defined such fhaf fhe Q facfor, defined as fhe signal efficiency divided by fhe square roof of fhe 
background efficiency, of a source wifh a specfral index of —2 is maximal by adjusfing fhe cuf on 
and fhe radius of fhe signal region 0 accordingly. A loose cuf sef is defined similarly but using 
an index of —3. An extraloose cut set designed for pulsar analyses is defined fo obfain maximal 
effeclive areas af low energies. Cuf values for and 0 are given in Tab. (1). 

All neural nefworks used in Ibis reconsfrucfion chain were frained for several observation 
conditions, differing in zenifh angle, azimufh angle and source offsef from fhe cenfre of fhe camera. 
In analyses of real dafa fhe MLP nefworks frained for simulafed positions adjacenf fo fhe real 
observafion position are used fo perform an inferpolafion. 

3. Performance of the Reconstruction Chain 

In fhe following fhe performance of fhe reconsfrucfion chain will be evaluated. All evenfs 
which have friggered CT 5 were analysed regardless of fhe frigger information of ofher telescopes 
in order fo demonsfrafe fhe performance over fhe entire coverable energy range. For example when 
performing splif observations wifh only CT 5 observing one source evenfs af all energies need fo 
be reconsfrucfed monoscopically. 

3.1 Direction Reconstruction Performance 

The direcfion reconsfrucfion performance can be sfudied by means of an angular resolution 
disfribufion as presented in Fig. (la) for fhree differenf zenifh angles. In Ibis figure fhe 68 % con- 
fainmenf radius Res of fhe disfribufion of reconsfrucfed direcfions around fhe direcfion of fhe simu¬ 
lafed 7 rays is shown as a function of fheir simulafed energy. The angular resolufion af low energies 
is ~ 0.3° and improves wifh increasing energy as fhe shower images become brighfer fo reach a 
minimum of ~ 0.1°, roughly comparable fo values achieved by fhe more advanced ImPACT anal¬ 
ysis, which uses a templale likelihood fif fo determine y-ray properfies [6], [7]. Towards highesf 
energies fhe angular resolufion worsens as fhe shower images become so large fhaf some images 
are cropped af fhe camera edges. In general fhis leads fo poorer even! reconstrucfion. 

3.2 Energy Reconstruction Performance 

The qualify of fhe energy reconsfrucfion can be quantified by fhe energy reconsfrucfion bias 
and resolufion. The bias is fhe mean and fhe resolufion is fhe RMS of fhe disfribufion of fhe relafive 
deviafion of fhe reconsfrucfed energies from fhe frue energy as a funcfion of fhe frue energy. Bofh 
quanfifies can, by definifion, only be calculated from Monfe Carlo dafa. 

The bias curves corresponding fo fhr'ee differenf zenifh angles are displayed in Fig. (lb). In ah 
cases fhe bias is smaller fhan 10% af energies up fo ~ 5TeV. Af low energies fhe bias is positive 
because only y-ray air showers whose recorded infensifies flucfuafed upwards can be defecfed. 

Energy resolution curves are shown in Fig. (Ic). Af fhe lowesf energies all shower images 
look approximately fhe same, fhus fhe reconsfrucfed energies are very similar leading fo a small 
resolution af lowesf energies. Af higher energies fhe resolufion sfays almosf consfanf around 20 - 
30%. Ah values are very similar fo fhe performance of fhe ImPACT analysis. 
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Figure 1: Instrument response functions corresponding to several zenith angles are shown in the 
above figures as functions of the true energy ETme- An azimuth angle of 0° and a wobble offset of 
0.5° was used to produce all curves. The displayed response functions are (a) the 68 % containment 
radius Rgg, (b) the energy bias, (c) the energy resolution and (d) the effective area. 


3.3 Effective Area 

The effective area is defined as the y-ray detection efficiency multiplied by the area corre¬ 
sponding to the maximum impact distance of the simulated air shower axes with respect to the 
telescope (500m at 20° zenith). Effective areas as a function of the true energy for three different 
zenith angles are shown in Eig. (Id). At the lowest energies the detection efficiency is very low. 
As the showers get brighter and more regular with increasing energy the efficiency increases until 
a maximum is reached. At the highest energies the efficiency begins to drop because the show¬ 
ers tend to get cropped at the camera edges leading to a worse classification performance. The 
maximum value of the effective areas increases with the zenith angle due to projection effects. 

3.4 Sensitivity 

The lowest point-source flux that can be detected with a significance of 5a within 50h of 
observations in an energy interval is called differential sensitivity. This quantity expressed in units 
of the differential Crab nebula flux is shown in Eig. (2) as a function of the reconstructed energy. 
The significances were estimated as N^/ ^NBkg, with being the number of excess events and 
Nskg being the exposure-corrected number of background events. The lowest detectable flux value 
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Figure 2: Differential sensitivities obtained by the MonoReco (black) and Model++ (red, [8]) anal¬ 
yses in units of the differential Crab Nebula flux as a function of the reconstructed energy Free- The 
curves correspond to a zenith angle of 20°, an azimuth angle of 180° and a wobble offset of 0.5°. 

per energy bin at a zenith angle of 20° corresponds to 4 % of the Crab Nebula flux. More advanced 
reconstruction chains like ImPACT or Model-i-i- [8], [9] reach values of ~ 3 %, as shown in Fig. (2). 

4. Analysis of Crab Nebula Data 

The reconstruction chain presented in the previous sections was applied to data from obser¬ 
vations of the Crab Nebula. In total 16 observation intervals lasting 28 min each were analysed, 
resulting in a live time of 7.2h. The observations were performed with a wobble offset of 0.5° 
from the camera centre at a mean zenith angle of 47°. Using the reflected background method [10] 
4800 excess events were measured in the ON region with respect to the number of events detected 
in OFF regions. This corresponds to a significance of 91a using Eq. (17) from Li & Ma [11]. 

The 0^ distributions shown on the left hand side of Fig. (3) visualise the distributions of the 
squared angular distances of the reconstructed shower directions from the centres of the corre¬ 
sponding ON and OFF regions. It can be seen that there is an excess towards the position of the 
Crab Nebula whereas there is none towards the test positions in the OFF regions. 

Significances can also be calculated from any position in the sky using the ring background 
method [10]. The resulting significance map is shown on the right hand side of Fig. (3). It can 
be seen that there are no artefacts visible in the field of view and that the source is point-like and 
centred on the position determined from radio measurements conducted with the VLBI [12]. 

At zenith angles around 45° the safe energy threshold is approximately 150GeV, compared to 
a thr'eshold of ~ 300GeV in case of H.E.S.S. I. The threshold is defined as the energy at which the 
effective area reaches 10% of its maximum value. This definition does not take the shape of the 
energy bias or energy resolution curves into account, thus the bias could, in principle, be arbitrarily 
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Figure 3: Left: Squared angular distance of the reconstructed directions to the test positions in the 
centre of the ON (black) and OFF (red) regions. Right: Sky map. 


large. Therefore a forward-folding technique was used to derive a differential energy spectrum. 
The resulting spectrum is shown in Fig. (4) together with the spectrum published by the MAGIC 
collaboration [13]. The spectral points are well-fit by a log-parabolic function. The fit parameters 
are given in the inset. The reference energy is 0.27 TeV. The spectrum is compatible with the result 
obtained by the MAGIC collaboration within systematic uncertainties. 

5. Summary and Outlook 

A new method of reconstructing the main properties of VHE y rays has been presented. The 
performance of the different parts of this reconstruction chain is competitive to other reconstruction 
concepts as discussed above. A big advantage of this chain is the speed with which events can be 
analysed. On average only 0.8 s of CPU time are required to analyse 5000 events, compared to 
approximately 20 s in the case of ImPACT. This makes the chain well-suited for an analysis of 
newly taken data on site in Namibia shortly after observations have been performed, enabling the 
H.E.S.S. collaboration to check for flux variabilify in fhe observed sources. Based on fhis, decisions 
can be made whefher fo confinue observing certain sources during fhe following nighfs. 

If will always be imporfanf fo have a monoscopic reconstrucfion chain available wifh which 
fhe entire energy range accessible to H.E.S.S. can be analysed, for example for split observation 
campaigns with only CT 5 observing a source. Still it is desirable to combine both monoscopic and 
stereoscopic reconstruction algorithms into one combined analysis, making it possible to decide on 
an event-by-event basis whether a monoscopic or a stereoscopic reconstruction algorithm should 
be used. Such an improvement will be provided in the near future. 
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